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Generalized Spin Glass Relaxation 
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Spin relaxation close to the glass temperature of CuMn and AuFe spin glasses is shown, by neutron 
spin echo, to follow a generalized exponential function which explicitly introduces hierarchically 
constrained dynamics and macroscopic interactions. The interaction parameter is directly related 
to the normalized Tsallis non-extensive entropy parameter, q, and exhibits universal scaling with 
reduced temperature. At the glass temperature q = 5/3 corresponding, within Tsallis' q-statistics, 
to a mathematically defined critical value for the onset of strong disorder and non-linear dynamics. 
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The dynamical correlations associated with the onset 
of a glassy state have been the focus of considerable in- 
terest and controversy over the last four decades, yet de- 
spite this attention such dynamics are still poorly un- 
derstood. Above the glass temperature, T g , the dynam- 
ics are generally characterized by a stretched exponential 
(Kohlrausch or KWW 0,01 ) relaxation, - exp[- (t/rf]. 
This form appears to be almost ubiquitous in nature, de- 
scribing phenomena as diverse as human dynamics 0, 
the conductivity close to the metal-insulator transition 
Q, the jamming transition frustrated magnets @, 
experimental [7] and theoretical [8[ spin glasses, as well 
as conventional glass forming liquids above the glass tem- 
perature 

However the stretched exponential form does not ad- 
equately describe the relaxation of either structural or 
magnetic glasses close to T g and below, where self- 
similarity in time (i.e. fractal behaviour in time) occurs. 
This is demonstrated in Monte Carlo calculations by 
Ogielski [8j], based upon a 3d ±J Ising spin glass model, 
which show that close to the spin glass temperature the 
time dependent spin autocorrelation function should take 
a phenomenological modified Kohlrausch form which in- 
corporates a power law dependence: 



q (t) = (Si (0) .Si (t)) oc t- x exp-^ 



(1) 



in which /3 increases from 1/3 at T 5 to 1 at ~ 4T gi x 
increases from below T g to 0.5 at high temperatures 
and the relaxation time, r, diverges at T g . This predic- 
tion is also largely supported by neutron spin echo (NSE) 
measurements [111 ]. 

The origins and implications of such non-exponential 
relaxation remain the subject of some debate, not least 
because it can be demonstrated that Kohlrausch-like re- 
laxation may arise from either a statistical distribution of 
independent (parallel) relaxation channels, or from more 
complex hierarchically constrained dynamics [l2[. For a 
detailed physical insight it is necessary to explore global 
models of glassy relaxation which attempt intrinsically to 



embody both distributed dynamics and the interactions 
that may lead to hierarchical relaxation. 

Such a model was introduced by Weron [13[ in an at- 
tempt to explain the apparently universal power law for 
dielectric relaxation. Weron's rigorous probabilistic ap- 
proach, based on the cluster model of Dissado and Hill 
[141 ], considers a hierarchical progression of relaxation 
which results in a continuously changing energy land- 
scape. The hierarchy is introduced through the formation 
of finite clusters which arise from interactions between 
relaxing dipoles, with each cluster being represented by 
an effective dipole related to its internal structure. The 
time taken for polarization fluctuations to reach equilib- 
rium is a random variable that for each relaxing dipole 
depends upon two other random variables, namely the 
waiting time and the dissipation rate. Through these 
variables Weron accounts for the effects of both inter- 
cluster and intra-cluster interactions with the character- 
istic timescale of any relaxing entity being restricted by 
the structural reorganization of the surrounding clusters 
and derives a generalized relaxation function 



<p(t) = 1 + k (t/r) 



-1/k 



(2) 



where /3 is associated with the fractal geometry of the 
system and k(> 0) is an effective interaction parameter 
related to the waiting time and providing a measure of 
the relative contribution of hierarchichal relaxation pro- 
cesses. Conveniently, the Weron power law reduces to 
the Kohlrausch form in the limit k — > 0, in which case 
< /3 < 1 has precisely the same meaning as before, with 
the limit (3 — > 1 implying simple Debye (exponential) re- 
laxation. 

Phenomenologically, the Weron model can be readily 
extended to the spin glass transition problem yet, despite 
its rigor and elegance, this approach has not previously 
been used in the analysis of spin relaxation. We have 
therefore investigated the applicability of this generalized 
relaxation function to spin glasses by analysing neutron 
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spin echo spectra collected from the archetypal metallic 
systems Cui-^Mn^ (x=0.1, 0.16, 0.35) and by revisit- 
ing our previously measured NSE spectra from Aui-^Fe^ 
(x=0.14)[Hj]. 

We have also explored the interconnections 
between the probabilistic Weron model and a more gen- 
eral approach based on non-extensive thermodynamics 
introduced by Tsallis [15] to describe highly disordered 
systems governed by Levy-stable distributions. 

Cui-^Mn^ samples with x=0.1, 0.16 and 0.35 were 
prepared by argon arc-melting and subsequent cold- 
rolling. The resulting disks were homogenised at 900° C 
before being quenched directly into water. The samples 
were found to have glass temperatures, T^, of 45K, 74K 
and 153K respectively. The Cui-^Mn^ NSE spectra were 
collected using the IN11C spectrometer [l6| at the Insti- 
tut Laue Langevin (Grenoble) at Q = 4nm _1 with an in- 
coming wavelength of 0.55nm. Our previously reported 
spectra from Auo.86Feo.14 (T^=41K) were obtained on 
IN15 at Q = 0.4nm _1 and a wavelength of 0.8nm [111 ] . 

All NSE spectra were independent of scattering vector, 
Q: Characteristic results are shown in figures [I] and O In 
figure [3] we show a direct comparison of least squares 
fits of the Kohlrausch stretched exponential, Ogielski 
(eqjl]) and Weron (eqj2|) functions to the relaxation of 
Auo.86Feo.14 at 45.7 K. Whilst both of the latter functions 
afford equally acceptable descriptions over the Fourier 
time range covered by the data, the Ogielski function 
must be modified to avoid an unphysical increase to val- 
ues greater than unity at short times [17]. In compari- 
son, the Weron function is able to describe all spectra, 
from the high temperature simple exponential limit to 
the complex power law decay found below T p , with con- 
sistent and physically meaningful parameters. 

The spectra of Auo.86Feo.14 give a temperature depen- 
dent /?, decreasing from 1 at high temperatures to ap- 
proximately 0.3 at the glass temperature, which implies 
a continuing evolution of the geometric and dynamical 
fractal character of the spin clusters. The CuMn spec- 
tra cover a smaller dynamic range, and equally valid fits 
can be obtained with a temperature independent /3 ~ 0.5 
or with (3 varying with temperature as found for AuFe. 
However, the most informative parameter is the interac- 
tion parameter, which diverges as the spin glass tem- 
perature is approached. We suggest that k provides an 
interesting and novel insight into the underlying ther- 
modynamics which drive the non-exponential spin glass 
relaxation, and closely links such relaxation to that ob- 
served in other disordered systems. 

A physical interpretation of the Weron function is pro- 
vided by Tsallis' generalisation of Boltzman-Gibbs ther- 
mostatistics for complex and multifractal systems [l5| . 
In this context Tsallis introduces the concept of non- 
extensive entropy, proposing that the total entropy of 
any self-organising, strongly interacting systems may be 
either greater or less than the sum of the entropies of the 
individual components of the system. 
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FIG. 1: Temperature dependence of the NSE spectra of (a) 
Cuo.9Mno.10 (b) Cuo.84Mno.i6 and (c) Cuo.65Mno.35 at Q = 
. In each case the solid lines represent fits of the Weron 
function, eqj2}, to the data. 

According to the Tsallis model, summing the entropy 
of two independent systems, A and £?, gives 

S q (A + B) _ S q (A) | S q (B) | S q (A)S q (B) 
ks ks ks ks ks 

where q is the so-called non-extensivity parameter. For 
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FIG. 2: Temperature dependence of the NSE spectra of 
Auo.86Feo.14 at Q = 0.4nra _1 . The solid lines are the as- 
sociated fits of the Weron function, eq(|2]). 
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FIG. 3: The NSE spectra of Auo.seFeo.M at 45. 7K with associ- 
ated least squares fits of the Kohlrausch, Ogielski and Weron 
relaxation functions described in the text. 



q > 1 the entropy is sub-extensive whilst for q < 1 it 
is super-extensive. Standard Boltzmann-Gibbs entropy 
(S = ks ^iPilnpi) is recovered as q —> 1. In general, the 
Tsallis total entropy may therefore be written as 



k B - 



(4) 



Brouers and Sotolongo-Costa [18[ demonstrated a di- 
rect relation between Tsallis' non-extensive entropy and 
Weron's probabilistic relaxation model, taking as their 
starting point the conventional weighted distribution of 
exponential relaxation processes 
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/ f(T)exp- t/T dT 
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(5) 



Here the relaxation time, r, scales with the characteris- 
tic volume v of the relaxing clusters according to r = v 1 ^ 
and < (3 < 1 is related to the fractal geometry of the 
system, as in the Weron analysis. By maximising the 
Tsallis non-extensive entropy they obtain a (normalised) 
cluster size distribution function 
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(6) 



where 1 < q < 2. This distribution function belongs 
to the family of Levy-stable distributions characterized 
by asymptotic power law tails [TH, [2(j and the non- 
extensivity parameter q quantifies the statistics of the 
system. An analysis of eq. [6] shows that the moments 
are defined only for l<q<3/2, whereas for q>3/2 all mo- 
ments, including the mean cluster size, diverge. How- 
ever, within the modified q- statistics introduced by Tsal- 
lis finite moments are obtained only for q<5/3 [21]. The 
critical value q=5/3 marks the transition to the limit of 
strong disorder, where most statistical weight is in the 
wings of the distribution and the macroscopic behavior 
is governed by a high number of highly improbable col- 
lective events. 

The relaxation function deduced from the cluster size 
distribution function takes a form identical to that pro- 
posed by Weron, subject to the simple substitutions 
k = (q-l)/(2-q) 



(7) 



The values of the interaction parameter, /c, obtained 
from fits of the Weron function of eq. [2] to the AuFe and 
CuMn NSE spectra can thus be related directly to the 
sub-extensivity parameter, q. 

As shown in figure [4j the resulting q- values for the all 
the spin glass alloys discussed here, together with those 
obtained from a re-analysis of the earlier published NSE 
measurements [22| on the more dilute Cuo.95Mno.05, ex- 
hibit the same scaling with reduced temperature, (T/T g ). 
At first sight this is surprising as the fractal character of 
dilute and concentrated CuMn spin glasses and of AuFe 
is expected to be quite different: Although an oscillatory 
RKKY magnetic exchange dominates the interactions in 
both systems, in AuFe the Fe atoms are known to cluster, 
whilst in CuMn the Mn atoms anticluster. Correspond- 
ingly ferromagnetic percolation is found at 15at%Fe in 
AuFe, but incipient antiferromagnetism is not evident in 
CuMn until 70at%Mn [23[. It is a remarkable prop- 
erty of phase transitions that critical phenomena do not 
depend on microscopic details but on the topology (i.e. 
interactions and dimensionality) of the systems. 

At high temperatures q is close to 1, indicating an es- 
sentially independent, parallel, relaxation of cluster mo- 
ments, q increases continuously with decreasing temper- 
ature and the first limit at q=3/2 is reached at ~ 1.2 
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FIG. 4: The Tsallis sub-extensivity parameter, q, obtained 
from the fits of eq(|2} to NSE spectra of the AuFe and CuMn 
samples, as function of reduced temperature. The values of 
q for the Cuo.95Mno.05 sample have been obtained by fitting 
the previously published data of ref [221 ]. 



T g . The most remarkable result, clearly demonstrated 
by figure 21 is that the spin glass transition temperature 
corresponds to the limit q ~ 5/3. If T^ marks the transi- 
tion to self-similarity, as suggested by MC simulations [8[ 
and experimental findings [ll| , then the relaxation at T^ 
should asymptotically follow the power law (t/r)~ y with 
y=/3x(2-q)/(q-l). Using /3 ~ 0.3 and q =5/3 we obtain 
an exponent y= 0.15, which is very close to y~0.13, the 
value found previously for AuFe and CuMn 5% [Tlj . 

The correspondence between T^ and a mathematically 
defined critical value for the extensivity parameter q, 
which defines a transition to complex non-linear dynam- 
ics, is new and has profound repercussions. It appears 
that spin glass systems evolve from a state characterised 
by conventional, extensive Boltzmann-Gibbs statistics 
at high temperatures to one of extreme sub-extensivity 
at T g . As q increases long-range correlations become 
more and more important and above q=3/2 the classical 
description breaks down, necessitating the introduction 
of the Tsallis q-statistics. The strong disorder limit is 
reached at q=5/3, where the complex dynamics are gov- 
erned by the power law tails of the Levy stable cluster 
size distribution function (eq. [6]). Large-amplitude ape- 
riodic fluctuations and nonlinearity dominate this limit, 
which can be identified with the (spin) glass transition 
and with the point, where self-similar relaxation sets- in. 



Significantly, this evolution appears to be universal, at 
least for the fundamentally different AuFe and CuMn 
spin glasses discussed here suggesting that spin glasses 
are very similar to many other complex disordered sys- 
tems such as financial markets, earthquakes, turbulence 
or jamming which are governed by by self-similarity and 
the underlining Levy stable distributions. 
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